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1.  Introduction 


This  technical  memorandum  studies  the  beamforming  of  large  aperture  vertical 
line  arrays  with  equally  spaced  sensors.  The  wavefield  sampled  by  a  vertical  array  is  the 
superposition  of  correlated  signal  arrivals  with  curved  wavefronts  in  a  spatially  colored 
ambient  noise  field.  The  interest  in  lower  frequencies  (10  to  100  Hz),  triggered  by  then- 
low  attenuation  in  the  water,  lead  to  the  development  of  low  frequency  acoustic  vertical 
line  arrays  such  as  the  MPL  digital  array  [Sotirinl988].  The  sampled  wavefield  is 
stror  ^iy  influenced  by  the  ocean  boundaries  since  the  physical  dimensions  of  such  arrays 
are  of  the  same  order  as  the  ocean  depth.  Also,  the  lower  frequency  sound,  which  has  a 
longer  acoustic  wavelength,  significantly  interacts  with  the  bottom.  One  objective  in 
developing  these  large  aperture  systems  is  to  improve  vertical  resolution.  An  increased 
resolution  allows  the  separation  of  the  multipath  arrivals  impinging  on  the  array,  thus 
making  possible  a  detailed  study  of  the  physics  of  sound  propagation.  Improved  resolu¬ 
tion  can  be  obtained  by  having  a  longer  aperture,  or  by  using  high  resolution  signal  pro¬ 
cessing  schemes  yielding  better  resolution,  or  both. 

Vertical  arrival  structure,  usually,  is  obtained  by  processing  the  data  sampled  by 
a  vertical  array  with  a  conventional  beamformer  using  plane  wave  replica  vectors.  Gen¬ 
erally,  the  assumptions  of  local  plane  wave  propagation  is  justified  by  the  moderate 
lengths  (a  few  hundred  meters)  of  the  apertures  deployed  in  the  water  column.  The  reso¬ 
lution  that  characterizes  this  conventional  processing  is  the  so-called  Rayleigh  resolution 
[Burdicl984]  given  by  0  (rad)  =  j-  where  X.  is  the  acoustic  wavelength  and  L  is  the  length 

of  the  aperture.  Increasing  the  aperture  length  results  in  an  improved  resolution  since 
9— >0  as  L  ->oo.  In  an  isovelocity  medium  with  sound  speed  of  1500  m/s,  the  1500  m 
long  array  described  in  [Sotirinl988],  has  a  Rayleigh  resolution  of  approximately  0.6°  at 
100  Hz  and  6°  at  10  Hz.  Even  for  such  a  long  array,  the  resolution  appears  limited  at  low 
frequency. 
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High  resolution  techniques  provide  a  substantial  improvement  of  resolution 
depending  on  the  encountered  signal-to-noise  ratio  SNR  .  High  resolution  techniques  can 
be  divided  in  two  classes,  the  adaptive  and  eigenvector  techniques.  In  the  following 
simulations,  the  adaptive  Minimum  Variance  Distortionless  Response  (MVDR)  beam- 
former  [Caponl969],  the  eigenvector  Multiple  Signal  Classification  (MUSIC)  algorithm 
[Schmidtl986]  and  the  Minimum  Norm  method  [Orfanidisl988]  are  considered. 

The  strong  potential  of  improved  resolution  with  these  high  resolution  tech¬ 
niques  also  results  in  poor  performance  as  soon  as  the  processed  wavefield  departs  from 
the  underlying  signal  model  on  which  they  are  based.  For  example,  the  high  resolution 
techniques  fail  to  report  correctly  correlated  arrivals  produced  by  multipath  propagator., 
which  is  common  in  a  real  oceanic  environment.  Spatial  smoothing  preprocessing  tech¬ 
niques  have  been  designed  to  limit  these  effects  at  the  cost  of  substantial  computational 
burden  and  of  reduced  effective  aperture  length  (thus,  a  reduction  in  resolution). 

The  complexity  of  sound  propagation  has  a  definite  impact  on  the  resolution  per¬ 
formance  of  the  beamformers.  Long  range  propagation  is  controlled  by  the  variations  of 
the  sound  speed  profile  with  depth  [Murphy  1987]  and  this  results  in  variations  of  the 
arrival  angle  with  depth  (wavefront  curvature).  This  curvature  can  cause  some  important 
smearing  and  bias  in  angle  in  the  arrival  structure  picture  produced  by  the  conventional 
processor  as  shown  in  [Tran  1989].  It  may  affect  the  high  resolution  methods  in  a  more 
critical  way  because  of  their  higher  sensitivity  to  mismatch.  Wavefront  Curvature  must 
be  taken  into  account  by  using  replica  vectors  that  properly  model  these  wavefronts. 

Coherent  processing  of  the  full  aperture  is  commonly  performed  using  the  con¬ 
ventional  beamformer,  which  is  a  very  computationally  efficient  processor  when  imple¬ 
mented  with  Fast  Fourier  Transforms  (FFTs).  It  is  natural  to  envision  high  resolution  full 
aperture  processing  of  a  long  and  well-filled  array,  although  several  aspects  of  this 
method  are  potentially  unpractical.  As  a  preliminary  step,  high  resolution  techniques 
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require  the  estimation  of  the  array  covariance  matrix.  In  practice,  the  array  covariance 
matrix  is  estimated  at  a  given  frequency  by  averaging  dyad  products  of  FFT  values  at 
each  sensor  over  successive  time  snapshots.  The  procedure  is  to  average  at  least  as  many 
dyads  as  sensors,  so  that  the  matrix  theoretically  is  of  full  rank.  Assuming  a  time 
snapshot  of  30  seconds  (0.033  Hz  FFT  bin  width)  and  an  array  of  120  elements,  the 
averaging  required  represents  an  hour’s  worth  of  data.  It  is  unlikely  that  the  wavefield 
will  be  stationary  for  that  long.  Furthermore,  the  high  resolution  techniques  involve 
computationally  expensive  matrix  operations  such  as  matrix  inversion  and  eigen  decom¬ 
position.  Coherent  processing  of  a  long  and  well  filled  array  requires  considerable  com¬ 
putational  resources  making  real-time  implementation  costly  and  difficult  [Speiserl985]. 

Beside  full  aperture  coherent  processing  with  high  resolution  techniques,  another 
strategy  with  high  resolution  methods  is  to  process  subapertures  extracted  from  the  verti¬ 
cal  array  and  to  incoherently  recombine  the  subaperture  angular  spectra.  Such  an 
approach  is  similar  conceptually  to  some  conventional  time  series  techniques  such  as  the 
Welch  periodogram  method  [Marplel987].  Even  though  using  smaller  apertures,  the 
subaperture  approach  should  yield  sufficient  resolution  due  to  the  incorporation  of  the 
high  resolution  methods,  and  should  enjoy  a  certain  level  of  robustness  due  to  the 
incoherent  averaging.  The  estimation  of  the  covariance  matrices  is  easier  with  their 
smaller  dimension  and  their  further  numerical  manipulation  facilitated.  In  addition  to  an 
angular  spectrum  estimate  across  the  entire  array  through  subarrays,  this  method  pro¬ 
duces  angular  spectra  corresponding  to  different  windows  in  the  water  column,  thus  pro¬ 
viding  information  related  to  multipath  variability  in  the  vertical. 

The  rest  of  this  technical  memorandum  does  an  intercomparison  of  the  two  dif¬ 
ferent  strategies  discussed  above: 

M )  a  full  aperture  coherent  processing  using  the  conventional  and  the  MVDR  beam- 

formers 
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(2)  a  subaperture  processing  using  the  high  resolution  techniques  (MVDR,  MUSIC 
and  the  Minimum  Norm  methods). 

The  analysis  is  based  on  processing  data  simulated  by  an  acoustic  modeling  program 
(based  on  a  normal  mode  decomposition)  that  outputs  the  complex  wavefield  at  range 
depth  pairs  of  interest.  The  simulated  wavefield  models  a  realistic  oceanic  environment 
with  correlated  and  curved  wavefront  arrivals. 
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2.  Creation  of  Synthetic  Data 

A  CW  tone  propagation  experiment  is  simulated  using  the  ATLAS  normal  mode 
model  [Gordonl984]  to  generate  the  complex  wavefield  which  outer  product  gives  the 
array  covariance  matrix  to  be  processed  by  the  various  methods.  Normal  mode  models 
provide  a  complete  solution  of  the  wave  equation  including  full  diffraction  effects 
[Boylesl984,  p.  199]  and  should  produce  a  realistic  wavefield  to  process  with  the  pro¬ 
posed  processing  methods  to  obtain  estimates  of  vertical  arrival  structure. 

2.1.  Description  of  the  Environment 

The  environment  is  characteristic  of  the  North-East  Pacific  and  a  deep  water 
situation.  It  is  range  independent  with  a  water  depth  of  4667  m.  The  sound  speed  profile 
with  depth  is  given  in  Table  2.1  and  plotted  in  Figure  2.1.  The  sound  speed  varies  in  the 
water  column  between  1489.1  m/s  and  1530  m/s.  The  sound  axis  is  500  m  deep.  The  crit¬ 
ical  depth  of  the  oceanic  waveguide  is  on  the  order  of  3313  m.  Thus,  the  sound  channel 
has  a  depth  excess  and  convergence  zone  propagation  may  take  place.  The  bottom  is 
characterized  by  a  frequency  invariant  bottom  loss  table  (bottom  loss  in  dB  as  a  function 
of  incidence  angle),  given  in  Table  2.2  and  plotted  in  Figure  2.2.  It  is  assumed  that  there 
is  no  surface  loss. 

A  5  m  deep  omnidirectional  noise  source  is  projecting  a  CW  tone  at  20  Hz  or 
100  Hz  with  a  source  level  of  200  dB  re  1  uPa  at  1  m.  The  sound  speed  at  the  source  is 
1509.46  m/s  so  that  its  reciprocal  depth  is  on  the  order  of  3275  m.  The  environment  is 
assumed  to  be  noise  free. 

When  the  frequency  analyzed  is  20  Hz,  the  receiving  vertical  array  is  assumed  to 
be  cut  for  30  Hz  with  an  interelement  spacing  of  25  m.  The  array  has  i23  identical 
equispaced  sensors  and  extends  across  most  of  the  water  column  from  100  m  to  3275  m. 
When  the  frequency  is  100  Hz,  the  receiving  vertical  array  is  assumed  cut  for  107  Hz 
with  an  interelement  spacing  of  7  m.  It  extends  from  300  m  to  1 189  m.  The  20  Hz  and 


6 


100  Hz  cases  differ  by  the  receiving  array  used.  The  geometry  in  each  case  is  summar¬ 
ized  in  Figure  2.3. 

2.2.  Analysis  and  Interpretation 

Before  processing  a  synthetic  wavefield  to  get  vertical  arrival  structure  informa¬ 
tion,  one  must  direct  one's  attention  toward  its  analysis  and  interpretation.  The  ray  trac¬ 
ing  approach  provides  an  easily  understandable  picture  of  the  wave  field  in  terms  of  ray 
diagrams,  and  guides  the  intuition  in  the  study  of  the  physics  of  the  simulation 
[Urickl983].  The  main  appeal  of  a  ray  tracing  program  is  the  direct  and  quick  availabil¬ 
ity  of  vertical  arrival  structure  information  across  the  array.  The  ray  tracing  program 
CONGRATS  in  the  Generic  Sonar  Model  (GSM)  will  be  used  in  the  following  study 
[Weinberg  1985], 

Ray  theory  is  widely  used  for  high  frequency,  deep  water  problems  [Jen- 
senl988].  As  rule  of  thumb,  one  may  use  in  practice  and  consider  the  ray  theory  solution 

to  be  fairly  accurate  if  the  frequency  /  is  such  that  /  >  where  c  is  the  reference 
sound  velocity  and  h  the  water  depth  [Comynl973].  This  is  not  an  absolute  criterion 
since  there  is  no  certainty  that  the  ray  solution  realistically  models  a  given  situation.  It  is 
well  known  that  ray  theory  is  a  high  frequency  approximation  to  the  solution  of  the  wave 
equation,  and  that  it  breaks  down  near  caustics  [Boviesl984],  For  this  reason,  rav  theory 
has  been  generalized  using  a  multipath  expansion  in  the  Fourier- Bessel  integral  represen¬ 
tation  of  the  pressure  field  to  numerically  solve  the  reduced  wave  equation  under  the 
WKB  approximation  [Weinbergl975].  This  modified  ray  theory  treats  the  caustic  prob¬ 
lem.  The  CONGRATS  program  in  the  GSM  implements  such  modified  ray  theory  and  is 
used  in  the  following.  It  is  considered  to  do  quite  well  at  low  frequency  in  its  computa¬ 
tion  of  eigenrays  and  pressure  [Weinberg  1985.  p.  6-123].  Pressure  or  transmission  loss  is 
estimated  by  a  phased  summation  of  the  different  eigenray  contributions  [Weinberg  1985. 
p.  6-139].  The  GSM  model  is  run  with  the  same  environmental  inputs  as  the  ATLAS 
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model.  The  joctom  loss  table.  Table  2.2,  corresponds  to  the  MGS  bottom  loss  table  Pro¬ 
vince  iype  3  of  the  GSM.  The  rays  with  launch  angles  of  0°,  ±5°,  ±  10°  and  ±  15°  from 
the  source  are  plotted  in  Figure  2.4. 

2.3.  20  Hz  Simulation  Results 

The  source  at  5  m  depth  projects  a  20  Hz  CW  tone.  The  transmission  loss  variations  with 
range  are  computed  by  the  ATLAS  model  and  the  results  plotted  on  Figure  2.5.  Panels 
A,  B  and  C  in  Figure  2.5,  respectively,  correspond  to  receiver  depths  of  5  m.  50  m  and 
500  m.  The  transmission  loss  at  a  range  of  80  km  is  120  dB  near  the  surface  and  on  the 
order  of  100  dB  at  greater  depth.  The  same  variations  of  transmission  loss  wuh  range  are 
computed  with  the  GSM  model  and  plotted  in  Figure  2.6.  The  GSM  results  have  the 
same  behavior  as  the  ATLAS  ones,  although  they  show  a  well  defined  convergence  zone 
between  50  and  60  km  which  is  not  reported  by  the  ATLAS  model.  The  acoustic 
wavelength  at  20  Hz  is  of  the  order  of  75  m  so  that  the  source  can  be  considered  at  the 
surface.  Both  models  report  consistent  transmission  loss  values  at  80  km  from  the 
source. 

The  receiving  array  is  assumed  to  be  at  56  km  from  the  source  where  the  conver¬ 
gence  zone  should  begin  if  it  exists.  The  variations  of  the  transmission  loss  as  a  function 
of  depth  across  the  128  sensor  array  computed  by  the  ATLAS  and  GSM  models  are 
respectively  plotted  in  Figure  2.7  and  Figure  2.8.  Although  the  interference  patterns 
across  the  water  column  do  not  match  exactly,  one  notes  the  fair  agreement  between  the 
two  models  below  KXX)  m  where  the  transmission  loss  varies  between  120  and  100  dB. 
As  noted  before,  the  GSM  results  yield  above  1000  m  a  convergence  zone  which  is  not 
reported  by  the  ATLAS  model.  The  phase  variations  across  the  array,  computed  by  the 
ATLAS  model,  are  plotted  in  Figure  2.9. 

The  arrival  angles  of  the  eigenrays  (rays  propagating  from  a  source  to  a  receiver) 
between  -45"  and  45"  tire  given  in  Table  2.3  as  calculated  by  the  GSM  at  56  km  range 
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and  at  the  sound  axis.  Eigenray  with  positive  source  angle  are  downward  going  at  the 
source  and  eigenrays  with  negative  target  angle  are  upward  going  [Weinuergi985,  p.  7- 
37],  In  order  to  present  results  consistent  with  the  beamforming  results  for  an  array  num¬ 
bered  from  top  to  bottom,  the  signs  of  target  angles  are  flipped  when  plotted,  like  in  Fig¬ 
ure  2.10  where  eigenray  angles  are  plotted  as  a  function  of  depth  at  a  range  of  56  km 
from  the  source.  In  Figure  2.10,  negative  angle  of  arrival  corresponds  to  downgoing  rav 
at  the  array  and  positive  angle  of  arrival  to  upgoing  ray  at  the  array.  One  observes  eleven 
major  arrivals  between  -45°  and  45°  at  the  array,  their  arrival  angles  at  the  array  vary 
significantly  with  depth.  The  directionality  picture  is  summarized  by  Figure  2.11  at  the 
sound  axis  depth  (500  m). 

2.4.  100  Hz  Simulation  Results 

The  CW  tone  projected  now  by  the  5  m  deep  source  has  a  frequency  of  100  Hz. 
The  receiving  array  is  cut  for  107  Hz  and  extends  from  300  m  to  1 189  m.  The  variations 
of  transmission  loss  with  range  computed  by  the  ATLAS  model  are  plotted  in  Figure 
2. 1 2  for  receiver  depths  of  5  m  (Panel  A),  50  m  (Panel  B)  and  500  m  (Panel  C).  Thest 
variations  me  similar  to  those  of  the  GSM  plotted  in  Figure  2.13.  The  better  agreement 
between  the  two  models  is  a  result  of  the  higher  frequency.  The  GSM  indicates  a  conver¬ 
gence  zone  between  50  and  60  km,  clearly  defined  for  a  receiver  at  50  m.  Approximately 
at  the  same  range,  the  transmission  loss  computed  by  ATLAS  decreases  slightly.  The 
convergence  zone  is  not  as  well  defined  as  in  the  GSM  results  because  low  frequency 
sound  diffracts  in  the  shadow  zone  [Boyles  1984,  p.  199-201].  The  existence  of  a  conver¬ 
gence  zone  is  confirmed  by  a  contour  plot  of  the  90  dB  transmission  loss  level  across  the 
whole  water  column  between  40  and  80  km  from  the  source  (Figure  2.14). 

Like  in  the  20  Hz  case,  the  receiving  array  is  assumed  to  be  at  56  km  from  the 
source,  in  the  far  field,  and  at  the  beginning  of  the  convergence  zone.  The  variations  of 
transmission  loss  across  the  128  sensors  are  plotted  in  Figure  2.15  and  2.16  for  the 
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ATLAS  and  GSM  results,  respectively.  The  transmission  loss  varies  rapidly  in  depth, 
indicating  a  complex  interference  pattern.  It  is  for  both  cases  on  the  order  of  90  to  95  dB. 
In  addition,  the  phase  variations  across  the  array,  computed  by  the  ATLAS  model,  are 
plotted  in  Figure  2.17. 

The  arrival  angles  of  the  GSM  eigenrays  between  -45°  and  45°  are  given  in 
Table  2.4  ,  as  calculated  by  the  GSM  at  the  sound  axis  and  for  a  range  of  56  km.  As 
above,  eigenrays  are  plotted  as  a  function  of  depth  with  their  signs  flipped  in  Figure  2.18 
to  be  consistent  with  the  sign  convention  that  corresponds  to  beamforming  an  array  num¬ 
bered  from  top  to  bottom.  One  observe  twelve  major  arrivals  between  -  45°  and  45°  on 
the  array,  their  arrival  angle  at  the  array  does  not  vary  significantly  with  depth  like  in  the 
20  Hz  case  since  the  aperture  is  three  times  smaller.  The  directionality  picture  is  sum¬ 
marized  by  Figure  2.19  at  the  sound  axis  depth  (500  m). 


Table  2.1  1 

Water  Deoth  (m) 

Sound  Velocity  (m/sl 
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Figure  2. 1 :  Sound  speed  profile. 
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Figure  2.5:  Transmission  loss  versus  range  computed  by  the  ATLAS  model  at  20  Hz. 
A:  receiver  depth  is  5  m,  B:  receiver  depth  is  50  m,  C:  receiver  depth  is  500  m. 
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Figure  2.6:  Transmission  loss  versus  range  computed  by  the  GSM  at  20  Hz.  A:  receiver 
depth  is  5  m,  B:  receiver  depth  is  50  m,  C:  receiver  depth  is  500  m. 


Table  2.3:  GSM  eigenray  table  at  20  Hz,  56  km  range  and  at  the  sound  axis. 


GENERIC  SONAR  MODEL  VERSION  ANSI 


RANGE  -  56  KM 

FREQUENCY  -  0.020  KHZ 

SOURCE  DEPTH  -  0.5C000e*01  M 
TARGET  DEPTH  -  0.50000*^-03  M 
FREQUENCY  -  0.20000e*-02  HZ 


ACOUSTIC  EIGENRAYS 


PANGE 

T  T  vir 

SOURCE 

TAPGET 

LEVEL 

FHASE 

n?pF 

\TB  T7 

KM 

s 

ANGLE 

DEG 
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DEG 

DB 
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56.0000 
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56 . 0000 
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1 

1 
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X 
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7 

1*; 
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1 V 

17 
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1 

: 
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11.8418 

-110.125 

-  ?  '  •<  - 

1 

\ 
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56 . 0000 

37 . 8044 
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2 
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2 
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-121.459 

-l'l .  03  l 

1 

2  v 
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27 
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2 
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2 

27 
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2 

27 
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2 

2 

56 . 0000 
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3 
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3 
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t 
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4 
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Figure  2. 10:  GSM  eigenray  arrival  angles  as  a  function  of  depth  at  -0  Hz. 
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Figure  2.12:  Transmission  loss  versus  range  computed  by  the  ATLAS  model  at  100 
A:  receiver  depth  is  5  m,  B:  receiver  depth  is  50  m.  C:  receiver  depth  is  500  m. 
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Figure  2.13:  Transmission  loss  versus  range  computed  by  the  GSM  at  100  Hz.  A: 
receiver  depth  is  5  m,  B:  receiver  depth  is  50  m,  C:  receiver  depth  is  500  m. 


Table  2.4:  GSM  eigenray  table  at  100  Hz,  56  km  range  and  at  the  sound  axis. 


GENERIC  SONAR  MODEL 


RANGE  -  56  KM 

FREQUENCY  -  100  H2 
SOURCE  DEPTH  -  5.0  M 
TARGET  DEPTH  -  503.0  M 

ACOUSTIC  EIGENRAYS 


RANGE 

TIME 

SOURCE 

TARGET 

LEVEL 

PHASE 

NSRF 

N8TM 

KM 

S 
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DEG 
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DEG 

DB 

DEG 

V  - 
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56, 0000 
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IV 
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-11.5864 
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IV 
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1 

1 

56 . 0000 
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1 

IV 
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1 

X  »' 

56.0000 
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IV 

3  '  T 
X  » 

56 . 0000 
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IV 

iv 
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n 

1 

56 . 0000 
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2 

IV 
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2 

:v 

56 . 0000 
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1 

2 

56 . 0000 
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1 
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]_ 

2v 
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-123 . 519 

-18“  .  517 

:v 

-7 
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;V 

27 
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-i 
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IV 
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-> 
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Figure  2.18:  GSM  eigenray  arrival  angles  as  a  function  of  depth  at  100  Hz. 
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3.  Beamforming  Results 


3.1.  Processor  Implementation 

The  two  complex  wavefields  created  by  the  ATLAS  normal  mode  program  are 
now  processed  by  the  full  aperture  and  subaperture  approaches  described  in  Section  1. 
The  conventional  beamformer  is  implemented  with  a  rectangular  window  and  a  Kaiser- 
Bessel  window  with  a  parameter  of  1.5  (yielding  a  38  dB  side  lobe  rejection).  The  angu¬ 
lar  spectra  are  normalized  to  yield  power  since  one  is  interested  in  estimating  discrete 
arrivals.  The  conventional  beamformer  using  plane  wavefront  replica  vectors  and  a  rec¬ 
tangular  window  is  given  by 

P rtct  =  ~2  Ep H  R  Ep  (3.1) 

where  the  plane  wavefront  replica  vector  Ep  is  such  that  EPHEP  =M.  M  is  the  number  of 
sensors  in  the  array  (it  is  equal  to  128).  The  normalization  factor  in  Equation  3.1 

produces  a  power  spectrum  (and  not  a  power  spectral  density).  R  is  the  array  covariance 
matrix  given  by 

R  =  X  Xw  (3.2) 

where  X  is  the  complex  wavefield  produced  by  the  ATLAS  normal  mode  model  at  the 
signal  frequency  of  interest.  When  the  Kaiser-Bessel  window  is  used,  the  conventional 
processor  associated  to  the  planar  array  manifold  is  given  by 


4/w=^(WE,)"R(WE,) 

where  Ep  and  R  are  defined  as  in  Equation  3.1  and  \V  is  defined  by 


(3.3) 


w  = 


l 

rnr= r 

XT  ^ 

M  i  -0 


a  o  0 

0  a  i 


6  6 


0 

0 

On-  i 


(3.4) 


where  the  a,’s  are  the  Kaiser-Bessel  window  coefficients  weighting  the  array,  M  is  the 
number  of  sensors.  Once  again  the  power  spectrum  produced  by  Equation  3.3  is 
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normalized  for  signals. 

The  phase  at  the  i,h  sensor  corresponding  to  a  plane  wavefront  replica  vector,  Ep , 
is  given  by 


~Zo)  (15) 
where  the  frequency  of  the  signal  is  / ,  0O  is  the  arrival  angle  with  respect  to  the  horizon¬ 
tal,  c0  the  sound  speed  (which  is  assumed  constant  throughout  the  whole  medium),  za  a 
reference  depth  and  z,  the  depth  of  the  i'h  sensor.  In  a  realistic  oceanic  medium,  the 
sound  speed  varies  with  depth.  The  sound  speed  profile  results  in  ray  bending  or  wave- 
front  curvature.  Even  for  modest  aperture  length,  such  curvature  can  be  important,  espe¬ 
cially  for  arrivals  with  small  angle  with  respect  to  the  horizontal  [Tran  1989],  Using 
Snell’s  law,  a  ray-geometric  approximation  of  the  wavefront  phase  was  derived  in 
[Tranl989],  The  phase  at  the  i'h  sensor  at  a  depth  z,  is  given  by 


A  f  Sgn(Q(zo)) 
l  ~  ci&)- 


1  -  COS20(Zo) 


c(Zi) 


c{Zo) 


dz 


(3.6) 


where  Sgn  is  the  signum  function,  and  cos26(z0) 


c(z.) 


C(Zo) 


<  l  is  implicitly  assumed  for  all  z, . 


As  earlier,  z0  is  the  reference  depth.  The  curved  wavefronts  phase,  which  enters  in  the 
curved  wavefront  replica  vector  Ec ,  exists  only  if 


ea’o^cos-1-^-  (3.7; 

If  Equation  3.7  is  not  satisfied,  the  corresponding  ray  has  turned  over  and  the  complex 
exponential  is  replaced  by  zero  at  the  sensor  position  i  in  the  vector  Ef.  Then,  partial 
insonification  takes  place  and  the  corresponding  curved  wavefront  replica  vector  Ec  is 


00 


0  •••  0 


(3.8) 


where  the  number  of  non  zero  entries  in  Equation  3.8,  genetically  denoted  by  * ,  is  A/,  ,  the 
size  of  the  insonified  aperture  (A/,  <M).  The  number  of  null  entries  and  their  position  in 
the  steering  vector  E0  depends  on  the  sound  speed  profile  and  the  particular  array 
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geometry  (i.e.  the  relative  position  of  the  array  with  respect  to  the  sound  speed  profile). 

The  conventional  beamformer  using  curved  wavefront  icplica  vector  is  given  by 


Wife?  E'"RE'  a9) 

where  the  insonified  window  size  A/,  depends  on  a  particular  curved  wavefront  steering 
vector  E,  or  look  direction,  E"  Ec  =  A/,  (EC).  A/;(EC)  is  generally  equal  to  M ,  except  for  low 
angle.  The  curved  wavefront  conventional  processor  using  a  Kaiser-Bessel  window  is 
similarly  given  by 


The  plane  and  curved  wavefront  replica  vectors  are  referenced  to  the  sound  axis 
so  that  all  the  results  obtained  by  the  various  processors,  operating  on  the  full  aperture  or 
on  subapertures,  are  consistent  with  each  other.  256  replica  vectors  between  -  90^  and 
+  90°  are  used,  the  quantization  in  angle  is  on  the  order  of  0.7° .  The  array  elements  are 
numbered  from  top  to  bottom  so  that  negative  angles  of  arrival  correspond  to  downgoing 
sound  or  uplooking  beams  and  positive  angles  of  arrival  correspond  to  upgoing  sound  or 
downlooking  beams.  This  sign  convention  is  consistent  with  the  results  presented  in 
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Section  2. 

In  addition  to  the  full  aperture  plane  and  curved  wavefront  conventional  beam¬ 
forming,  the  MVDR  processor  is  used  to  beamform  the  full  aperture  as  well  as  subaper¬ 
tures  extracted  from  the  128-sensor  array.  The  MVDR  beamformer,  in  the  following,  is 
given  by 


Psn'DR  e,«  m),-1  Ef  (3'13) 

where  F(R),  is  the  spatially  smoothed  covariance  matrix  on  s  -sensor  subarrays  across  a 
M  sensor  array  (the  spatial  smoothing  subsegment  length  is  s ).  When  the  full  aperture  is 
processed,  M  is  equal  to  128  and  the  spatial  smoothing  subsegment  length  s  is  64  sensors. 
When  subaperture  processing  is  performed,  M  is  chosen  equal  to  32  and  the  spatial 
smoothing  subsegment  length  j  equals  20.  The  curved  wavefront  steering  vectors  associ¬ 
ated  to  each  M -sensor  array  (each  one  has  a  particular  geometry  with  respect  to  the  sound 
speed  profile)  are  passed  through  the  spatial  smoothing  pre-processing  operation  from 
which  the  "best”  replica  vectors  are  obtained  by  an  eigenvalue-eigenvector  decomposi¬ 
tion  as  explained  in  [Tranl989]  ;  hence  the  name  of  eigensteering  vector  Ef.  The 
eigensteering  vector  E, ,  of  dimension  (j,1  ),  is  normalized  to  unit  norm.  Er  11 E,  =  1.  In  the 
case  of  partial  insonification,  the  eigensteering  vectors  corresponding  to  low  angle  of 
arrival  with  respect  to  the  horizontal  are  given  by 


E, 


(3.14) 


where  the  number  of  non  zero  entries  in  Equation  3.13  generically  denoted  by  *  is  r,  (E, ). 
the  size  of  the  insonified  equivalent  aperture  resulting  from  spatial  smoothing  with  length 
s .  As  noted  before,  r,(E,),  generally,  is  equal  to  s  except  for  low  angle  of  look.  The  nor¬ 
malization  by  s,  ( E, )  in  Equation  3.12  is  performed  to  yield  a  power  spectrum  and  not  a 
power  spectral  density.  F(R),  is  stabilized  before  inversion  by  adding  a  small  fraction  of 
noise  corresponding  to  10-4  of  the  average  power  across  the  s  sensor  equivalent  smoothed 
array  (i.e.  lO-VtEtR),),  where  ir  denotes  the  trace  operation).  This  allows  for  a  condition 


36 


number  (ratio  of  the  smallest  to  the  largest  eigenvalue)  on  the  order  of  10-6. 

The  MUSIC  and  Minimum  Norm  eigenvector  methods  are  used  to  do  subaper¬ 
ture  processing,  like  the  MVDR  beamformer.  The  MVDR,  MUSIC  and  Minimum  Norm 
methods  will  operate  on  32-element  subapertures  using  the  corresponding  eigensteering 
vectors  for  spatial  smoothing  with  subsegment  length  of  20  sensors.  Subaperture  process¬ 
ing  is  performed  on  48  subarrays  sliding  from  the  top  to  the  bottom  of  the  full  aperture 
array  by  two  elements  (i.e.  the  first  array  extends  from  sensor  1  to  32,  the  second  subar¬ 
ray  extends  from  sensor  3  to  34  and  so  on).  The  angular  spectra  are  computed  on  each 
subaperture  and  then  incoherently  averaged  to  produce  a  full  aperture  composite  angular 
spectrum. 

The  subaperture  processing  is  also  performed  on  four  non-overlapping  32  sensor 
subarrays,  using  the  corresponding  eigensteering  vectors  as  described  above.  The  compo¬ 
site  angular  spectrum,  then,  is  the  results  of  averaging  on  those  four  subarrays.  Such  pro¬ 
cessing  appears  more  realistic,  being  more  practical  in  its  computational  requirements. 

The  MUSIC  algorithm  is  given  by 

P"USIC  =  TTE7T  e,"  rL,  e,  t3-15) 

where  s;(E,)  and  E,  are  defined  as  before.  is  the  noise  only  covariance  matrix.  It  is 
obtained  by  choosing  a  number  of  arrivals  or  signals  n,ignol  and  doing  an  eigenvalue- 
eigenvector  decomposition  of  the  spatially  smoothed  covariance  matrix  F<R>:.  R,,,..,  is 
given  by 


R„ 


l"~  V,  V/' 


(3.16) 


S  ^  signal  1  *  1 

where  the  V,  "s  are  the  unit  norm  noise  eigenvectors  (V,  V,w=  l).  The  orthonormal  eigen¬ 
vectors  V,  of  F(R),  corresponds  to  the  j  ~n„sml  lowest  eigenvalues.  To  ensure  numerical 
stability  in  the  computation  of  the  quadratic  form  in  Equation  (3.15).  the  noise  only 
covariance  matrix  is  stabilized  by  adding  a  small  fraction  of  noise  to  main  diagonal,  i.e. 
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lO~7tr(RMije),  where  tr  denotes  the  trace  operation.  This  results  in  clipping  the  peaks  SNR 
theoretically  infinite  to  70  dB. 

The  Minimum  Norm  method  is  given  in  [Orfanidisl988,  p.  363] 


p si, norm  =  eTWe,  0A1) 

where  Ee  is  a  (s,l)  unit  norm  eigensteering  vector  as  defined  earlier,  and  d  a  ($.1)  vector 
given  by 


d  = 


.2  V.V« 


Uo 


(3.18) 


where  U<f  =  [l  0  0  •  ■  •  0],  and  V,  like  in  MUSIC  is  the  /'*  noise  eigenvector  of  the  spatially 
smoothed  covariance  matrix  F(R),.  To  ensure  numerical  stability  of  the  computation  of 
the  quadratic  form  in  Equation  3.7,  the  (s,  .y)  matrix  ddw  is  stabilized  by  adding  a.  small 

fraction  of  white  noise  to  the  main  diagonal,  lO-7— Idd'iL,  where  ir  denotes  the  trace 

operation.  As  in  the  case  of  the  MUSIC  algorithm,  this  results  in  clipping  theoretically 
infinite  peaks  SNR  to  70  dB. 
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3.2.  20  Hz  Simulation  Results 

The  results  of  conventional  processing  on  the  20  Hz  synthetic  data  using  rec¬ 
tangular  and  Kaiser-Bessel  windows,  and  the  plane  wave  array  manifold  are  plotted  in 
Figures  3.1.  These  angular  spectra  correspond  to  the  output  of  a  classic  FFT  beamformer 
commonly  performed  on  vertical  line  arrays.  The  situation  is  of  high  signal-to-noise  ratio 
since  there  is  no  ambient  noise.  The  highest  power  arrivals  lay  between  -  30°  and  +  30° . 
The  use  of  a  tapering  window  appears  necessary  to  protect  the  estimate  from  sidelobe 
leakage,  the  high  angle  power  level  drops  off  20  dB  when  one  switches  from  the  rec¬ 
tangular  window  to  the  Kaiser-Bessel  window. 

The  use  of  curved  wavefront  replica  vectors  with  the  conventional  beamformer 
gives  results  similar  to  the  case  of  plane  wavefront  replica  vectors.  The  angular  spectra 
for  the  curved  wavefront  array  manifold  are  plotted  in  Figures  3.2  for  rectangular  and 
Kaiser-Bessel  windows.  The  arrivals  are  reported  at  slightly  different  angles.  The  con¬ 
ventional  beamformer  using  curved  wavefront  replica  vectors  reports  high  power  levels 
on  the  order  of  90  dB  at  near  horizontal  angles.  This  is  due  to  the  partial  insonification  of 
the  aperture  which  results  in  higher  sidelobes  and  lower  resolution.  The  insonified  aper¬ 
ture  in  the  128  element  array  is  plotted  as  a  function  of  arrival  angle  in  Figure  3.3.  Near 
horizontal  arrivals  insonify  a  very  limited  part  of  the  aperture,  thus  explaining  the  high 
power  levels  reported  at  low  angle  by  the  curved  wavefront  processor.  Horizontal  arrivals 
in  this  case  do  not  physically  exist  since  they  insonify  only  the  sensor  number  13  as 
shown  in  Figure  3.3. 

The  results  of  the  conventional  beamformer  using  plane  and  curved  wavefront 
replica  vectors  are  overlaid  with  the  GSM  eigenrays  in  Figure  3.4.  The  GSM  multipath 
pattern  at  the  sound  axis  is  similar  to  the  angular  spectra  obtained  from  the  ATLAS  data 
with  all  arrival  angles  steeper  than  ±10".  The  peak  powers  indicated  by  the  GSM  are 
slightly  higher  than  the  ones  in  the  ATLAS  angular  spectra.  This  is  consistent  with  the 
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acoustic  modeling  performed  on  Section  2  which  shows  that  the  convergence  zone 
reported  by  the  GSM  is  much  better  defined  than  the  ATLAS  one.  One  can  note  a  slight 
increase  in  resolution  when  curved  wavefront  replica  vectors  are  used,  especially  for  low 
arrival  angles.  There  is  also  a  better  agreement  with  the  GSM  eigenray  arrival  angles  for 
the  two  arrivals  near  -15°  and  20°  for  the  curved  processor.  Furthermore,  the  curved 
wavefront  processor  yields  arrival  peak  to  in-between  peak  levels  larger  than  the  plane 
wavefront  processor  for  physically  upgoing  arrivals  with  angles  up  to  +40°  and  also 
downgoing  arrivals  with  arrival  angles  larger  than  -30°.  The  upgoing  arrival  angular  esti¬ 
mates  produced  by  the  curved  processor  are  also  closer  to  the  GSM  eigenray  arrival 
angle  than  the  ones  produced  by  the  plane  wave  processor.  This  is  also  true  for  the 
downgoing  arrival  close  to  -15°.  The  planar  replica  vectors  produce  slightly  lower 
arrival  angles,  especially  for  upgoing  sound  (positive  angles).  The  conventional  proces¬ 
sor  produces  similar  results  for  steep  arrival  angles  with  either  type  of  replica  vectors.  In 
this  case,  it  is  known  than  the  phase  relationships  entering  in  the  plane  and  curved  wave- 
front  replica  vectors  are  almost  equivalent  [Tran  1989]. 

The  output  of  the  MVDR  processor  operating  on  the  full  aperture  with  a  spatial 
smoothing  subsegment  length  of  64  sensors  is  plotted  in  the  top  panel  of  Figure  3.5.  A 
large  amount  of  smoothing  is  necessary  in  order  to  decorrelate  over  a  dozen  multipath 
arrivals  which  can  be  close  together,  as  indicated  by  the  GSM  eigenray  angles.  There¬ 
fore,  the  MVDR  estimate  may  still  suffers  large  signal  cancellation  due  to  arrival  correla¬ 
tion.  The  peak  power  of  the  arrivals  is  up  to  20  dB  lower  than  the  ones  produced  by  the 
conventional  processor.  Loss  can  also  be  due  to  incorrect  wavefront  modeling,  even  if 
eigensteering  vectors  are  used,  and  unsufficient  sampling  in  angle.  The  MVDR  beam- 
former  is  known  to  be  very  sensitive  to  mismatch,  especially  at  high  signal-to-noise  ratio. 
The  near  horizontal  angles  (at  0°)  exhibits  an  increase  in  power  which  is  due,  like  in  the 
case  of  the  conventional  curved  wavefront  processor,  to  the  partial  insonification  of  the 
spatially  smoothed  equivalent  aperture.  The  eigensteering  vectors  which  corresponds  to 
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this  equivalent  aperture  have,  at  low  arrival  angles,  zero  entries  much  like  the  actual  phy¬ 
sical  aperture,  as  shown  in  Figure  3.6.  The  results  of  the  MVDR  full  aperture  processing 
are  overlaid  with  the  curved  wavefront  conventional  processing  results  and  the  GSM 
eigenrays  in  the  bottom  panel  of  Figure  3.5.  Although  the  equivalent  aperture  is  reduced 
by  half  after  spatial  smoothing,  the  MVDR  beamformer,  using  the  eigensteering  vectors, 
yields  a  much  larger  resolution  than  the  conventional  beamformer. 

The  results  of  the  32  sensor  subaperture  processing  by  the  MVDR  method,  after 
spatial  smoothing  with  a  subsegment  length  of  20  and  using  the  eigensteering  vectors,  are 
presented  in  Figures  3.7,  3.8  and  3.9.  Figure  3.7  is  the  waterfall  plot  of  the  angular  spec¬ 
tra  corresponding  to  each  subarray  and  referenced  to  the  sound  axis.  Increasing  subarray 
number  corresponds  to  increasing  depth.  This  plot  shows  that  the  number  of  peaks,  as 
well  as  their  power,  are  not  correctly  reported.  Because  of  the  short  aperture,  only  limited 
spatial  smoothing  can  be  performed.  This  leads  to  an  unsuffcient  decorrelation  of  the 
multipaths  and  signal  cancellation  takes  place.  The  flat  area  at  low  angles  in  the  waterfall 
plot  corresponds  to  an  angular  region  where  no  arrivals  or  sound  can  physically  exist 
according  to  the  ray  geometric  approach  selected  here.  It  is  the  result  of  the  partial 
insonification  phenomena  discussed  earlier.  In  Figure  3.8  are  plotted  the  traces  with 
depth  of  the  local  maxima  of  these  spectra,  each  angular  spectrum  is  associated  to  the 
depth  of  the  corresponding  subarray  center.  A  peak  is  detected  and  indicated  in  Figure 
3.8  if  it  is  a  local  maximum  in  the  angular  spectrum  (null  derivative  and  negative  second 
derivative).  This  plot  shows  the  great  variability  of  the  detected  arrival  angles  with  depth 
(even  though  referenced  to  the  sound  axis).  It  is  largest  in  the  middle  and  the  upper  half 
of  the  128  sensor  aperture.  As  a  result  of  this  great  variability,  the  incoherent  averaging 
of  the  48  angular  spectra  (all  referenced  to  the  sound  axis)  does  not  produce  good  results, 
as  shown  in  Figure  3.9.  The  arrivals  have  peak  powers  which  are  greatly  lower  than  the 
conventional  beamformer  peak  powers.  In  addition  the  arrival  structure  is  significantly 
different  from  the  one  obtained  by  full  aperture  coherent  processing  with  either  the 
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conventional  beamformer  or  the  MVDR  processor. 

The  same  processing  is  performed  by  the  MUSIC  algorithm.  The  number  of  sig¬ 
nals  is  12.  The  MVDR  and  MUSIC  methods  are  known  to  yield  similar  results  at  very 
high  signal  to  noise  ratios  [Nickel  1988].  Since  the  situation  studied  is  ambient  noise 
free,  it  is  the  case  here  and  the  MUSIC  angular  "spectra"  waterfall  plot  of  Figure  3.10  is 
qualitatively  similar  to  the  MVDR  angular  spectra  waterfall  plot  of  Figure  3.7.  The 
traces  with  depth  of  the  local  maxima  of  these  spectra,  plotted  in  Figure  3.11,  show  the 
same  variability  as  in  the  case  of  the  MVDR  beamformer.  The  incoherent  average  of  the 
MUSIC  "spectra",  plotted  in  Figure  3.12,  is  almost  identical  to  the  MVDR  results  in  Fig¬ 
ure  3.9,  with  an  imperceptible  improved  resolution.  This  can  be  observed  easily  by  look¬ 
ing  at  Figure  3.13  on  which  are  overlaid  the  average  of  the  48  subaperture  spectra  for  the 
MVDR  beamformer  and  the  MUSIC  algorithm,  and  the  GSM  eigenray  angles.  The  spec¬ 
tra  are  normalized  to  their  respective  minimum  value.  One  observes  a  good  agreement  in 
angle  estimation  between  the  MUSIC  averaged  "spectrum"  and  the  GSM  results  for  low 
arrival  angles.  However,  the  angular  estimate  differ  for  steeper  angle  of  arrival  which 
also  corresponds  to  lower  powers. 

The  Minimum  Norm  method  is  used  with  the  same  spatial  smoothing  subseg¬ 
ment  length  of  20  and  a  number  of  signals  set  to  12.  The  waterfall  plot  of  the  Minimum 
Norm  "spectra"  in  Figure  3.14  shows  a  substantial  increase  in  resolution  compared  to  the 
MUSIC  results,  with  more  defined  peaks.  The  traces  with  depth  of  the  local  maxima  of 
these  "spectra"  are  plotted  in  Figure  3.15.  The  variability  is  even  larger  than  with 
MUSIC.  As  a  result,  the  incoherent  averaging  of  the  48  Minimum  Norm  "spectra"  pro¬ 
duces  a  complicated  picture  with  a  high  density  of  peaks  in  Figure  3.16.  The  MUSIC 
and  Minimum  Norm  averaged  "spectra"  are  normalized  by  their  respective  minimum 
value  and  overlaid  with  the  GSM  eigenray  arrival  angles  in  Figure  3.17.  There  is  a  good 
agreement  in  angle  estimation  at  low  angles  of  arrival.  The  Minimum  Norm  estimate  has 
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line  pairs  in  place  of  single  arrival.  This  appears  similar  to  the  line  splitting  phenomenon 
common  for  high  resolution  spectral  estimation  methods  such  the  auto-regressive  model 
based  methods  [Marplel987]. 

The  averaged  angular  spectra  based  on  four  nonoverlapping  subarrays  of  32  sen¬ 
sors  (with  a  spatial  smoothing  with  subsegment  length  of  20)  are  plotted  in  Figures  3.18, 
3.19  and  3.20  for  the  MVDR,  MUSIC  and  Minimum  Norm  methods,  respectively.  The 
composite  angular  spectrum  is  the  average  of  the  spectra  calculated  over  four  non¬ 
overlapping  32-sensor  subarray.  Averaging  over  four  subarrays  in  some  cases  produces 
better  results  than  averaging  over  48  subarrays,  as  done  earlier,  because  of  the  great  vari¬ 
ability  of  the  angular  spectra  in  the  vertical.  The  improvements  are  characterized  by  a 
better  resolution  of  the  peaks,  and  are  especially  obvious  for  the  Minimum  Norm  results 
where  the  multipath  arrival  structure  is  fairly  close  to  the  eigenray  structure  produced  by 
the  GSM,  The  MUSIC  and  Minimum  Norm  results  are  overlaid  with  the  GSM  eigenray 
angles  in  Figure  3.21.  One  notes  for  both  MUSIC  and  the  Minimum  Norm  the  qualitative 
improvement. 
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3.3.  100  Hz  Simulation  Results 

The  100  Hz  case  is  characterized  by  an  aperture  length  smaller  than  in  the  20  Hz 
case,  over  which  sound  speed  does  not  vary  as  much.  Therefore,  the  effects  of  the  wave- 
front  curvature  should  have  less  impact.  The  angular  spectra  produced  by  the  conven¬ 
tional  beamformer,  using  planar  replica  vectors,  with  a  rectangular  window  and  a 
Kaiser-Bessel  window,  are  plotted  in  Figure  3.22.  The  angular  spectra  produced  by  the 
conventional  processor  using  curved  wavefront  replica  vectors  are  plotted  in  Figure  3.23. 
As  in  the  20  Hz  case,  a  tapering  window  is  necessary  to  minimize  sidelobe  leakage.  The 
curved  wavefront  conventional  processor  produces  high  power  in  the  horizontal  because 
of  the  partial  insonification  of  the  aperture,  which  results  in  large  sidelobe  leakage 
[Tran  1989].  The  insonified  aperture  is  plotted  in  Figure  3.24  which  shows  that  the  0° 
arrival  insonifies  only  three  sensors. 

The  results  of  the  conventional  processing  with  plane  and  curved  wavefronts  are 
overlaid  with  the  GSM  eigenrays  in  Figure  3.25.  As  expected,  the  conventional  proces¬ 
sor,  using  curved  and  plane  wavefront  replica  vectors,  are  identical,  except  for  the  spuri¬ 
ous  horizontal  arrival  of  the  curved  wavefront  processor.  They  match  closely  the  mul¬ 
tipath  arrival  structure  reported  by  the  GSM.  The  correspondence  between  the  ATLAS 
angular  spectrum  estimate  and  the  GSM  eigenrays  is  still  excellent  when  looking  at  the 
absolute  power  levels  of  the  arrivals. 

The  output  of  the  MVDR  processor  operating  on  the  full  aperture  array,  after 
spatial  smoothing  with  a  subsegment  length  of  64  sensors,  is  plotted  in  the  top  panel  of 
Figure  3.26.  The  MVDR  angular  spectrum  is  overlaid  with  the  results  of  the  conventional 
processor  using  curved  wavefront  replica  vectors  and  the  GSM  eigenrays  in  the  bottom 
panel  of  Figure  3.26.  As  in  the  20  Hz  case,  the  MVDR  processor  on  the  full  aperture 
method,  using  eigensteering  vectors,  produces  a  better  resolution  than  the  conventional 
processor,  even  with  a  reduction  by  half  of  the  aperture  length  due  to  spatial  smoothing. 
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The  peak  power  levels  in  this  angular  spectrum  are  up  to  20  dB  lower  than  the  ones  of 
the  conventional  processing.  There  may  still  be  some  strong  signal  cancellation  in  addi¬ 
tion  to  mismatch  loss  due  to  imperfect  wavefront  modeling  and  unsufficient  sampling  in 
angle.  The  slight  increase  of  power,  at  low  arrival  angles,  corresponds  to  the  partial 
insonification  of  the  spatially  smoothed  equivalent  aperture.  The  msonified  part  of  this 
equivalent  aperture,  based  on  the  eigensteering  vectors  is  plotted  in  Figure  3.27. 

The  results  of  the  32  sensor  subaperture  processing  by  the  MVDR  method,  after 
spatial  smoothing  with  subsegment  length  of  20  and  using  the  eigensteering  vectors,  are 
presented  in  Figures  3.28,  3.29  and  3.30.  Figure  3.28  is  a  waterfall  plot  of  the  angular 
spectrum  corresponding  to  each  subarray.  As  in  the  20  Hz  case,  these  spectra  are  refer¬ 
enced  to  the  sound  axis.  The  number  of  peaks  in  Figure  3.28  is  roughly  half  of  what  is 
detected  in  the  full  aperture  coherent  processing.  This  is  due  to  the  limited  spatial 
smoothing  performed  on  the  short  subaperture.  Here,  spatial  smoothing  does  not  fully 
decorrelate  the  multipaths.  Thus,  signal  cancellation  takes  place  and  results  in  lower 
power  levels.  Another  cause  of  peak  power  reduction  is  mismatch  due  to  imperfect  wave- 
front  modeling.  The  peaks  levels  are  up  to  20  dB  below  the  conventional  processor  peak 
levels.  The  traces  with  depth  of  the  local  maxima  of  these  spectra  are  plotted  in  Figure 
3.29.  The  high  power  level  and  low  angle  arrivals  are  stable,  while  the  lower  level,  high 
angle  arrivals  exhibit  a  large  variability  both  in  power  and  angle.  The  composite  angular 
spectrum,  calculated  by  incoherently  averaging  the  48  angular  spectra,  is  plotted  in  Fig¬ 
ure  3.30.  The  strong  and  low  angle  arrivals  detected  by  this  composite  estimator  have 
power  levels  on  the  same  order  as  the  full  aperture  MVDR  processor. 

The  subaperture  processing,  performed  by  the  MVDR  processor,  is  now  repeated 
using  the  MUSIC  algorithm,  with  a  number  of  signals  set  to  12.  The  waterfall  plot  of  the 
MUSIC  "spectrum"  in  Figure  3.31  shows  some  resolution  improvement.  The  traces  with 
depth  of  the  local  maxima  of  these  "spectra",  plotted  in  Figure  3.32,  still  indicate  a  large 
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variability  for  steeper  arrivals.  The  limited  spatial  smoothing  does  not  allow  decorrela¬ 
tion  of  all  the  multipath  arrivals.  13  spatial  averages  are  performed  to  transform  the  array 
covariance  matrix.  Based  on  the  output  of  full  aperture  processors,  at  least  15  multipath 
arrivals  impinge  on  the  array.  Incoherent  averaging  of  the  MUSIC  "spectra"  in  Figure 
3.33  produces  a  very  good  estimate  of  the  arrival  structure  for  arrival  angles  between 
-  35°  and  +  35° ,  even  if  it  does  not  have  as  high  signal  to  floor  levels  as  in  the  case  of  the 
full  aperture  processors.  For  steeper  arrival  angles,  which  have  also  lower  power  levels, 
the  MUSIC  algorithm  fails  to  report  the  correct  arrival  angles. 

The  subaperture  composite  spectra  for  the  MVDR  and  MUSIC  are  normalized 
by  their  respective  minimum  value  and  overlaid  with  the  GSM  eigenray  angles  in  Figure 
3.34.  Both  processors  detect  the  highest  level  arrivals  with  no  or  little  angle  bias  at  low 
arrival  angle.  The  bias  in  angle  is  larger  for  steeper  and  lower  level  arrivals.  As  suggested 
earlier,  this  is  due  to  insufficient  decorrelation  achieved  by  limited  spatial  smoothing. 

The  Minimum  Norm  method  now  is  used  to  perform  the  subaperture  processing 
performed  by  the  MVDR  and  MUSIC  techniques,  with  a  20  sensor  subsegment  length 
spatial  smoothing  and  a  number  of  signals  set  to  12.  The  waterfall  plot  of  the  Minimum 
Norm  angular  "spectrum"  with  depth,  in  Figure  3.35,  shows  a  substantial  increase  in 
resolution,  with  peaks  more  defined.  The  traces  with  depth  of  the  local  maxima  of  these 
"spectra"  are  plotted  in  Figure  3.36  and  exhibit  a  larger  variability  for  higher  arrival 
angles.  On  the  other  hand  the  lower  angle  arrivals  are  very  stable  with  depth.  As  a  result, 
the  average  of  the  "spectra"  estimated  on  the  48  subarrays,  plotted  in  Figure  3.37,  has 
strong  and  well  defined  peaks  at  low  angles,  and  more  choppy  noise-like  structure  at 
higher  angles.  This  is  due  to  unsufficient  decorrelation  achieved  by  the  spatial  smoothing 
with  a  subsegment  length  of  20,  which  only  allows  the  resolution  of  the  strongest 
arrivals. 

The  averaged  Minimum  Norm  and  MUSIC  "spectra"  are  normalized  by  their 
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respective  minimum  value  and  overlaid  with  the  GSM  eigenray  angles  in  Figure  3.38. 
We  observe  a  very  good  agreement  between  the  GSM  eigenray  angles  and  the  arrival 
angle  estimates  obtained  from  the  MUSIC  and  the  Minimum  Norm  techniques.  The 
Minimum  Norm  method  allows  an  easier  detection  of  arrivals  compared  to  the  MUSIC 
algorithm,  yielding  larger  signal-to-noise  ratios. 

The  subaperture  processing  on  four  32  sensor  nonoverlapping  subarrays  is  now 
presented  in  Figures  3.39,  3.40  and  3.41  for  the  MVDR,  MUSIC  and  Minimum  Norm 
techniques  respectively.  Both  the  MVDR  composite  spectrum  and  the  MUSIC  compo¬ 
site  "spectrum"  based  on  the  4  subarrays  are  similar  to  the  composite  spectra  based  on  the 
48  subarrays.  It  appears  in  Figures  3.39  and  3.40  that  doing  48  averages  or  only  four 
averages  has  no  or  little  effects  on  the  detection  of  the  main  low  angle  arrivals.  On  the 
other  hand,  the  Minimum  Norm  composite  "spectrum",  based  on  only  4  non-overlapping 
subarrays,  allows  the  detection  of  additional  arrivals  which  impinge  on  the  array  with 
steeper  angles,  compared  to  the  MUSIC  algorithm  (Figure  3.42). 
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4.  Conclusions 

This  report  introduces  several  approaches  to  beamform  large  vertical  line  arrays. 
Full  aperture  beamforming  using  the  conventional  and  the  MVDR  processors  is  pro¬ 
posed.  Conventional  beamforming  using  plane  and  curved  wavefront  replica  vectors  is 
performed.  Full  aperture  MVDR  processing  is  accomplished  using  curved  wavefront 
replica  vectors  which  are  modified  to  minimize  mismatch  after  spatial  smoothing.  A 
subaperture  processing  is  described  and  implemented  using  three  high  resolution 
methods,  the  MVDR,  MUSIC  and  Minimum  Norm  techniques  with  curved  wavefront 
replica  vectors  modified  to  minimize  mismatch  after  spatial  smoothing. 

A  simulation  is  designed  to  assess  the  relative  performance  of  these  various  pro¬ 
cessing  schemes.  It  does  not  include  ambient  noise  and  thus  corresponds  to  infinitely 
high  signal-to-noise  ratios.  The  ATLAS  normal  mode  model  is  used  to  create  realistic 
environment  signals  to  be  processed.  Careful  modeling  allows  a  good  correspondence  to 
be  achieved  between  the  ATLAS  beamforming  results  and  the  GSM  eigenrays,  with  only 
marginal  differences  in  angle  of  arrival  and  power  levels.  Twc  cases  are  treated.  The  20 
Hz  case  is  of  a  very  large  array  operating  at  low  frequency,  where  almo  the  whole 
water  column  is  covered.  The  100  Hz  case  corresponds  to  mid-low  frequencies  and 
arrays  currently  available  such  as  the  MPL  digital  array  [Sotirinl988]. 

The  use  of  curved  wavefront  replica  vectors  with  the  conventional  beamformer 
appears  to  be  desirable  in  the  20  Hz  case.  The  use  of  plane  wavefront  replica  vector  can 
result  in  arrival  angle  bias  in  arrival  angle,  especially  for  low  angle  of  arrival  with  respect 
to  the  horizontal  where  curvature  is  the  largest  [Tran  1989].  The  use  of  curved  wavefront 
also  yield  in  some  cases  better  peak  resolution.  In  the  100  Hz  case  where  the  array  is  900 
m  long,  the  simulation  shows  that  there  is  no  advantage  of  using  curved  wavefront 
replica  vectors  with  the  conventional  processor.  In  this  case,  the  output  of  the  conven¬ 
tional  beamformer  is  almost  the  same  whether  planar  or  curved  wavefront  replica  vectors 
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are  used. 

Because  the  simulation  does  not  include  any  ambient  noise  model,  the  signal-to- 
noise  ratios  are  infinitely  high  and  loss  due  to  mismatch  can  become  a  limiting  factor  for 
the  performance  of  the  high  resolution  methods.  The  use  of  curved  wavefront  replica 
vectors,  modified  to  accommodate  the  spatial  smoothing  pre-processing  (which  is 
required  in  order  to  decorrelate  the  multipaths),  is  necessary  in  both  the  20  Hz  and  the 
100  Hz  cases.  The  MVDR  processor  is  used  to  process  coherently  the  full  aperture  and 
produces  better  resolution  than  conventional  beamforming.  On  the  other  hand,  loss  due 
to  mismatch  and  signal  cancellation  due  to  multipath  arrivals  still  occur  and  result  in 
lower  power  levels  than  in  conventional  beamforming.  The  stabilization  noise  added  to 
the  smoothed  covariance  matrix  in  order  to  invert  it,  sets  a  noise  floor  in  the  angular 
spectra. 

The  subaperture  processing,  which  requires  less  integration  time  to  estimate  the 
array  covariance  matrix,  does  not  produce  as  good  signal-to-noise-flocr  levels  as  the  full 
aperture  processing.  Based  on  the  simulation,  the  high  resolution  eigenvector  based 
methods.  MUSIC  and  the  Minimum  Norm,  have  fairly  good  resolution  capabilities,  even 
if  the  parameters  used  in  this  work  are  such  that  the  amount  of  spatial  smoothing  is 
potentially  insufficient  to  decorrelate  all  the  multipath  arrivals. 

The  simulation  brings  up  the  problem  of  selecting  the  location  and  the  number  of 
subarrays  to  process  in  order  to  obtain  the  averaged  angular  spectrum  of  arrival.  The 
vertical  arrival  structure  appears  highly  variable  depending  on  the  amount  of  incoherent 
averaging  performed  (e.g.  when  48  subarrays  or  4  nonoverlapping  subarrays  me  pro¬ 
cessed  by  the  MUSIC  and  the  Minimum  Norm  methods).  Increasing  the  number  of 
subarrays  does  not  necessarily  enhance  the  angular  spectrum.  Too  much  incoherent 
averaging  may  result  in  smearing  the  peaks  in  the  angular  spectrum. 

This  simulation  shows  that  the  two  processing  strategies  envisioned  are  viable 
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and  that  the  high  resolution  methods  using  a  ray  geometric  model  of  curvature,  based  on 
the  sound  speed  profile  at  the  array,  produce  good  results.  By  using  the  high  resolution 
methods,  one  enhances  resolution  in  the  angular  spectra.  The  cost  to  pay  in  this  simula- 
uon  is  a  reduction  of  signal-to-noise-floor  level  essentially  because  of  the  absence  of 
ambient  noise. 
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Figure  3.1:  20  Hz  output  of  the  conventional  beamformer  using  plane  wavefront  re¬ 
plica  vectors.  Data  are  weighted  in  the  top  panel  by  the  rectangular  window,  and  in 
the  bottom  panel  by  a  Kaiser-Bessel  (o  =  1.5)  window. 
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Figure  3.2:  20  Hz  output  of  the  conventional  beamformer  using  curved  wavefront 
replica  vectors.  Data  are  weighted  in  the  top  panel  by  the  rectangular  window,  and 
in  the  bottom  panel  by  a  Kaiser-Bessel  (n  -  1.5)  window. 
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Figure  3.4:  20  Hz  output  of  the  conventional  'ueamformer  with  a  Kaiser-Bessel  win¬ 
dow  with  n  parameter  of  1.5  using  curved  wavefront  replica  vectors  (solid  line), 
plane  wavefront  replica  vectors  (dotted  line)  overlaid  with  the  GSM  eigenrays  at  the 
sound  axis. 


Figure  3.7:  Waterfall  plot  of  the  32  sensor  subarray  MVDR  angular  spectra  at  20  Hz 
after  spatial  smoothing  with  subsegment  length  of  20  sensors  and  using  the 
eigensteering  vectors. 
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Figure  3.10:  Waterfall  plot  of  the  32  sensor  subarray  MUSIC  "spectra"  at  20  Hz  after 
spatial  smoothing  with  a  subsegment  length  of  20  sensors,  using  the  eigensteering 
vectors,  and  with  a  number  of  signals  set  to  12. 
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Figure  3.14:  Waterfall  plot  of  the  32  sensor  subarray  Minimum  Norm  "spectra"  at  20 
Hz  after  spatial  smoothing  with  subsegment  length  of  20  sensors  using  the 
eigensteerlng  vectors  and  with  a  number  of  signals  set  to  12. 
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Figure  3.19:  MUSIC  "spectrum"  (average  on  4  non-overlapping  32  sensor  subarray 
spectra  at  20  Hz).  Spatial  smoothing  is  performed  with  subsegment  length  of  20 
sensors.  The  number  of  signals  is  set  to  12. 
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Figure  3.20:  Minimum  Norm  "spectrum"  (average  on  4  non-overlapping  32  sensor 
subarray  spectra  at  20  Hz).  Spatial  smoothing  is  performed  with  subsegment  length 
of  20  sensors.  The  number  of  signals  is  set  to  12. 
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Figure  3.21:  MUSIC  (solid  line)  and  Minimum  Norm  (dotted  line)  subaperture  "spec¬ 
tra"  (averaged  on  4  non-overlapping  32  sensor  subarrays)  overlaid  with  the  GSM 
elgenray  angles  at  the  sound  axis  and  at  20  Hz. 
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Figure  3.22:  100  Hz  output  of  the  conventional  beamformer  using  plane  wavefront 
replica  vectors.  Data  are  weighted  by  the  rectangular  window  In  the  top  panel  and 
by  the  Kaiser-Bessel  («  =  1.5)  window  in  the  bottom  panel. 
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Figure  3.23:  100  Hz  output  of  the  conventional  beamformer  using  curved  wavefront 
replica  vectors.  Data  are  weighted  by  the  rectangular  window  in  the  top  panel  and 
by  the  Kaiser-Bessel  window  or  =  1.5  in  the  bottom  panel. 
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Figure  3.25:  Angular  spectra  at  100  Hz  produced  by  the  conventional  beamformer 
with  a  Kaiser-Bessel  window  with  n  parameter  of  1.5  using  curved  wavefront  replica 
vectors  (solid  line),  plane  wavefront  replica  vectors  (dotted  line)  overlaid  with  the 
GSM  elgenrays  at  the  sound  axis. 
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Figure  3.26:  Top  panel:  100  Hz  output  of  the  full  aperture  MVDR  beamformer.  Bot¬ 
tom  panel:  overlaid  MVDR  (solid  line)  and  conventional  (dotted  line)  spectra  and 
GSM  eigenrays  at  the  sound  axis. 


Figure  3.27:  Insonlflcatlon  of  the  64-sensor  spatially  smo 
by  the  elgensteering  replica  vectors  between  -90”  and  90*. 
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Figure  3.28:  Waterfall  plot  of  the  32  sensor  subarray  MVDR  angular  spectra  at 
100  Hz  after  spatial  smoothing  with  subsegment  length  of  20  sensors  and  using  the 
eigensteering  vectors. 
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Figurp  3.30:  Composite  MVDR  angular  spectrum  at  100  Mz. 
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Figure  3.31:  Waterfall  plot  of  the  32  sensor  subarray  MUSIC  "spectra"  at  100  Hz 
after  spatial  smoothing  with  a  subsegment  length  of  20  sensors,  using  the 
eigensteerlng  vectors,  and  with  a  number  of  signals  set  to  12. 
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Figure  3.35:  Waterfall  plot  of  the  32  sensor  subarray  Minimum  Norm  "spectra"  at 
100  Hz  after  spatial  smoothing  with  subsegment  length  of  20  sensors  using  the 
eigensteering  vectors  and  with  a  number  of  signals  set  to  12. 
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Figure  3.39.  MVDR  subaperture  angular  spectrum  (averaged  on  4  non-overlapping 
32  sensor  subarray  spectra  at  100  Hz).  Spatial  smoothing  is  performed  with  a  sub- 
segment  length  of  20  sensors. 
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Figure  3.41:  Minimum  Norm  subaperture  "spectrum"  (averages  on  4  non¬ 
overlapping  32  sensor  subarray  spectra  at  100  Hz).  Spatial  smoothing  is  performed 
with  subsegment  length  of  20  sensors.  The  number  of  signals  is  set  to  12. 


92 


BIBLIOGRAPHY 


Boyles,  A.  C.,  Acoustic  Waveguides,  Applications  to  Oceanic  Sciences,  John  Wile7, 
N.Y.,  1984. 

Capon  J.,  “High  Resolution  Frequency-Wavenumber  Spectrum  Analysis,”  Proc.  IEEE , 
’  vol.  57,  pp.  1408-1418,  1969. 

Corny n,  J.  J.,  GRASS:  A  Digital  Computer  Ray  Tracing  and  Transmission  Loss  Predic¬ 
tion  System,  1,  NRL  Report  7621,  21  December  1973. 

Gordon,  D.  F.  and  H.  P.  Bucker,  Arctic  Acoustic  Propagation  Model  with  Ice  Scattering, 
NOSC  Technical  Report  985,  30  September  1984. 

Jensen,  F.  B.,  “Wave  Theory  Modeling:  A  Convenient  Approach  to  CW  and  Pulse  Pro¬ 
pagation  Modeling  in  Low-Frequency  Acoustics,”  IEEE  ,J.  Ocean.  Eng.,  vol.  13, 
no.  4,  pp.  186-197,  October  1988. 

Marple,  S.  L.,  Digital  Spectral  Analysis  with  Applications,  Prentice-Hall,  N.J.,  1987. 

Murphy,  D.  A.  and  D.  R.  DelBazo,  “Multipath  Processing  for  Improved  Detection  on  a 
Long  Vertical  Array,”  in  Progress  in  Underwater  Acoustics,  pp.  765-772,  Ple¬ 
num  Press,  N.Y.,  1987. 

Nickel,  U.,  “Algebraic  Formulation  of  the  Kumaresan-Tufts  Superresolution  Method 
Showing  Relation  to  ME  and  MUSIC  Methods,”  IEE  Proc.  Part  F,  vol.  135.  no. 
l,pp.  7-10,  1988. 


Orfanidis,  S.  J.,  Optimum  Signal  Processing,  An  Introduction,  Macmillan.  New  York, 
1988. 

Schmidt,  R.  O.,  “Multiple  Emitter  Location  and  Signal  Parameter  Estimation.”  Trans. 
Ant.  Prop.,  vol.  34,  no.  3,  pp.  276-280,  1986  . 

Sotirin,  B.  J.  and  J.  A.  Hildebrand,  “Large  Aperture  Digital  Acoustic  Array."  IEEE  J. 
Oceanic  Eng.,  vol.  13,  no.  4,  pp.  271-281,  October  1988. 

Speiser,  J.  M.,  “Progress  in  Eigenvector  Beamforming,”  SPIE  Proceedings,  vol.  564, 
San  Diego,  CA,  August  22-23. 

Tran,  J-M.  and  W.  S.  Hodgkiss,  “High  Resolution  Beamforming  on  Vertical  Arrays  in  a 
Realistic  Oceanic  Environment,”  MPL  Technical  Memorandum  408,  MPL-U- 
03/89,  Marine  Physical  Laboratory,  Scripps  Institution  of  Oceanography,  San 
Diego,  CA,  January  1989  . 

Urick,  R.  J.,  Principles  of  Underwater  Sound,  McGraw-Hill.  N.Y.,  1983. 

Weinberg,  H.,  “A  Continuous-Gradient  Curved  Fitting  Technique  for  Acoustic-Rav 
Analysis,” ./.  Acotist.  Soc.  Am.,  vol.  50,  no.  3  Part  2,  pp.  975-984.  March  1971 . 


93 


Weinberg,  H.,  Generic  Sonar  Model,  NUSC  Technical  Document  597 1-D,  6  June  1985. 


